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(57) Abstract 

A system for locating objects in space, such as medical instniments (16) within the body of a patient (P), based upon transmission of 
magnetic fields from coils (10) in a fixed frame of reference to sensors (20. 22. 24) on the object or vice-vetsa. The cunent supplied to the 
coils is adjusted to assure that the sensors receive fields within a preselected range of magnitudes regardless of ^ location of the object 
in space. Hiis assures that the smsor operates within its optimum range, and pennits use of compact transmitters and s^sors. 
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DEBCRTPTTQiy 

MAGNETIC LOCATION SYSTEM 
WITH ADAP TIVE PEEDBACK GOMTOQT, 

TECHNIGAI, FIEI,D 

The present invention relates to systems for 
determining the location and/ or orientation of objects 
in space by detecting magnetic fields. 
BACKGROUND ART 

Various systems have been proposed for 
detecting the position and/or orientation of an object 
using magnetic or electromagnetic fields. These systems 
typically employ field transmitters, such as 
electromagnet coils, disposed at known locations in a 
fixed reference frame and a sensor, such as a coil or 
other transducer mounted to the object to be located. 
Each transmitter projects a field varying in space in a 
fixed frame of reference. The pattern of variation in 
space for each transmitter is different than the pattern 
for each other transmitter. For example the 

transmitters may be identical to one another but 
disposed at different locations or in different 
orientations. The field patterns of the transmitters 
are thus displaced or rotated relative to one another 
and relative to the fixed frame of reference. The sensor 
on the object detects the parameters of the field 
prevailing at the location of the object as, for 
example, the magnitude and/ or direction of the field at 
the object or the magnitudes of individual components of 
the field at the object in one or more preselected 
directions.. The transmitters may be actuated in a 
predetermined sequence so that at any time only one 
transmitter is active and therefore the field prevailing 
at the object is only the field contributed by one 
transmitter, plus a background field due to the Earth's 
magnetic field and other environmental sources. 
Alternatively, the transmitters can be driven at 
different frequencies so that components of the signal 
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from the sensor varying at different frequencies 
represent contributions to the field at the object fro. 

different transmitters. Based «««« \ ^ 

. eased upon the detected 

5 trir °' individual 

tran=„xtt.rs, and th. taown pattern or variation ot the 
fxeld irom each transmitter, a computer syst«. 
calculates the position and orientation of the ej^eor 
and hence the position of the object hearin, the .^Z] 
xn the £«ed frame of reference of the transmitters, 

10 a variant of this system, the ohject to he loLa^ 
M^alL - transmitters, whereas^ 

or or\.Tta: "T" 

or orientations in the fi«a frame of reference. The 

15 f«m s" T'" « Object is deduced 

» from signals representing the parameters of the field 
prevailing at the various sensors. 

in united states Patents 4,849,692, «,64J 786- 
4,710,708; 4,613,866 and 4 945 5n6 = 1 •-142,786, 
20 4,945,305. Systems according 

» to this general design can be used to provide a three! 
dimensional spatial input capability for a computer, 
another system of this nature is disclos^ In 
international Patent publication W094/04938. i„ the 

2S medi r r"°"' "^"^^ ^ ^ «>e « 

r = include a sensor 

mounted on the tip of an endoscope, so that the location 
and/or orientation of the endoscope tip can b2 

det«min.d While the sensor is disposed insi^ «r boS 
Of a medical patient. This allows the physician to 
30 monitor the endoscopic procedure without resorti:, to ! 
fluoroscopy or other techniques using ionizing radiation 
11 T . ^"'^"^^^ on- embodim«.t (pages 

in Which T r""*'°" - arran^LL 

35 ,? " transmitting coils are actuat«l 

35 simultaneously to "steer» the direction of the result!^ 

u!Ll.t, " Characterizes this arrangement as 

undesirable. other systems for locating medical 
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instruiaent:s such as endoscopes and catiheters based upon 
transmitted fields are disclosed in United States Patent 
No, 5,042,486; 5,099,845; 5,211,165; 5,251,635; 
5,253,647; 5,255,680; 5,265,610 and 5,391,199» 

Systems of this nature typically have used 
fields having strength which varies as the third power, 
or a higher power, of distance from the transmitter • 
There is, accordingly, a very large variation in the 
strength of each field from location to location* In a 
medical device locating system, adapted to detect the 
location of the sensor anywhere within a sensing volume 
having dimensions on the order of 0,5 to l.o m, the 
field strength may vary through many orders of magnitude 
from one end of the sensing volume to another. When the 
sensor and the object to be located happen to be near a 
particular transmitter, they will receive an extremely 
strong field, whereas when the sensor and object are at 
the end of the sensing voliune remote from the particular 
transmitter, they will receive only a very weak field 
from that transmitter. This requires a sensor with a 
very substantial dynamic range, capable of accurately 
monitoring very strong and very weak fields. That, in 
turn, poses stringent demands on the design of the 
sensor and may require a larger sensor. Moreover, the 
signal to noise ratio is poor when monitoring a very 
weak field and the accuracy of the system is accordingly 
less than optimum in these conditions. Both of these 
considerations are particularly important in systems for 
medical applications, where the sensor must be small to 
fit inside a medical device such as a catheter or 
endoscope. Although such systems can provide useful 
results, further improvement would be desirable. The 
same problems arise in systems using a coil or other 
transmitter on the object to be located and plural 
sensors in the fixed frame of reference. Here again, a 
particular sensor will be exposed either to a very 
strong field or a very weak field depending on whether 
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5 ^„v„ "'^^'2,479 discloses a solution to these 

^: .7; T''" ~^^">'» accoram, 

the .479 application include electro»a,„ets which are 

fields. Each field has at least one component with a 
non-zero magnitude that is either constant T 
nearl. linear with respect to distamra pa^i^i:: 
dxrection within a sensing volume such P^"^«^lar 
provides less variation L flTiH arrangement 
^, , a^Tiaiiion xn fxeld magnitude throughout 

whLThe WelfT " ""^""^ arran,e„e„r " 

wnich the field strength varies as the thlra ^ v., ^ 

» power Of distance. Thus, with the ^asi^^^ta^ ^U?^ 
the difference between .i„i^ Txl^u^^I^d 
a given sensing volume is substantially s«.ller. 
suh»t^tially alleviates the problems discussed 
However, it is not always convenient to use th. 

» ^rticular coil arrangements which provide :uch ;.aS! 
linearly varying fields. 

=.n„,v. V ""^ "^"^ °* invention provides 

waratus for determining position including Weld 
^ generating means for producing a plurality of magnet" 
fields. The field generating means is arrang^ to 
provide such fields having parameters, such as fieU 

variation. The pattern of variation for each field is 
d fferent than «,e pattern of variation for eac^ othe^ 
field Por example, wh«:e th. field generating means 
includes a plurality of transmitting colls, th! 
are disposed at different locations and/or different 
35 orxentations relative to <-ho a^- • airrerent 
accordina t« 4.K sensing volume. Apparatus 

at T!ast "^'"^ °' invention also includes 

at least one sensor adapted to detect one or more 
parameters o. the .ield prevailing at the sensT. 7Z 
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the sensor is at an unknown location within the sensing 
volume, and to provide one or more sensor signals 
representing such detected parameters. For example, the 
sensor may include a sensor body and a plurality of 
5 component sensors disposed on the sensor body, each 
component sensor being operative to measure the 
magnitude of a magnetic field component in a preselected 
local direction relative to the sensor body to provide a 
component sensor signal representing that particular 
component of the field prevailing at the sensor. The 
apparatus further includes calculation means for 
calculating the location of the sensor, the orientation 
of the sensor or both based upon the sensor signals and 
upon the unknown patterns of variation of the fields. 
The apparatus also includes feedback control means for 
adjusting the field generating means to alter the known 
pattern of variation of at least one of the fields 
responsive to the sensor signals, to the calculated 
location of the sensor or both so as to maintain the 
detected parameters of the altered field at the sensor 
within a preselected range. For example, where the 
field generating means includes a plurality of 
transmitters disposed adjacent the sensing volume, the 
feedback control means can be arranged to increase the 
25 strength of the field transmitted by a particular 
transmitter when the sensor is remote from such 
transmitter, and to decrease the strength of the field 
emitted by that transmitter when the sensor is close to 
the transmitter. Thus, the parameters of the field 
30 detected at the location of the sensor will always lie 
within a relatively narrow range of values. Because the 
transmitted field is altered in a known manner, the 
position and/or orientation of the sensor can still be 
calculated, simply by using the altered field as the 
35 basis for calculation. For example, where the feedback 
control means increases the current flowing to an 
electromagnetic coil, the parameters used in the 
calculation to represent the strength of the field 
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emitted from that particular r.«4 1 

accordingly. ''"^^ adjusted 

Apparatus according to this aspect of the 

present invention can maintain the detected fie 1^ 

Of the field at such location, within a relatively 
narrow range even where the field strength varies with I 
power of distance greater than the first power aJ^or 
example, where the field strength varies with the th!rl 
power Of distance from the coil. The system thtr^^or! 
can provide good signal to noise ratio, even wittTv::; 
small sensor. Moreover, there is no need for a^ 
particular coil configuration to produce a linear :r 
quasx-lxnear field. Por example, in a medical 

applxcatxon, a plurality of individual coils may he 
dxsclosed at various locations beneath a patient 
receiving bed, or on one side of the bed. P^^^«nt 
Although apparatus according to this aspect of 

additional advantages when used with magnetoresistive 
sensors or other sensors which can lose accuracy w^In 
exposed to very high magnetic fields. Because the flel^ 
xs maintained within a relatively narrow range, the 
sensors are not exposed to fields which would Lpair 

thrf'drr* ' '""^^^^ "^^'^"^ invention, 

the feedback control system is actuated to maintain th^ 
parameters of the field at each sensor within ^e 

30 ! maintains the field parameters at the 

30 sensor at substantially constant values. Here aga!^? 
the transmitted fields are altered in a 3.nown maLr 
and the disposition of the sensor can still bl 

A further aspect of the present invention 
aZrr /"""^'"^ trin^itt:" 
also including a plurality of sensors disposed adjacent 



20 



25 



35 



wo 96/41119 



PCTAJS96/08411 



a sensing volvune. The apparatus according to this 
aspect of the invention includes calculation means for 
determining a calculated location of the transmitter 
based upon signals from the sensors and upon the known 
5 pattern of variation of the field at the transmitter. 
Here again, feedback control means are provided for 
adjusting the field generating means to alter the known 
pattern of variation of field sent by the at least one 
transmitter in response to the sensor signal, to the 

10 calculated location of the transmitter or both so as to 
maintain the parameters of the field from the at least 
one transmitter prevailing at each sensor within a 
preselected range. This arrangement is similar to the 
arrangement discussed above, except that the transmitter 

15 or transmitters are disposed on the object to be located 
whereas the sensors are located in the fixed frame of 
reference of the sensing volxme. Here again, when the 
transmitter is close to the sensor, or where the 
orientations of the transmitter and sensor favor strong 

20 coupling between the transmitted field and the sensor, 
the strength of the field is reduced. If the 

transmitter is relatively far from the sensor, or in a 
unfavorable orientation, the control means increases the 
strength of the transmitted field. This arrangement 

25 provides advantages similar to those discussed above. 

A further aspect of the invention provides 
apparatus for sensing the disposition of an object — its 
position, orientation or both — in a frame of reference. 
Apparatus according to this aspect of the invention 

30 includes transmitter means including at least one 
transmitter for providing a field and sensor means 
including at least one sensor for sensing one or more 
parameters of a field prevailing at each such sensor and 
providing one or more sensor signals indicative of such 

35 parameters. The transmitter means and said sensor means 
cooperatively define a plurality of transmitter-sensor 
pairs, each including one transmitter and one sensor as 
elements of the pair, one element of each such pair 
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being disposed on the oblect anri 

«^ wjjjecT: ana the other elein«an4- 

5 different . ^^"^^^ Pair is disposed at a 

Within a preselected range. 

provide ::::z t'^i^l^ "-^^^n 

« orientation w.erein the rr?t\r„, "r'nene^X 

r^s^rr^eirnr,?.^^^^^^ 

Of tho ensor, or from the calculated location 

::ra:::.rr;;u' rrrr- ^ '-'-^"^ - 

» Within a p.e.electVrai^e ^ 

Other objects, features and advantages of tho 
IZIT.T"^'"" "'^^^ appare„r.r:l 

f^ti lei v"°" °' ""'^"-^ embodiments 

« a::l:e ' --l™-^- -^^^ ao=o.^„nn, 

d^nlof ^^^"".^ ^ ^ diagramiaatic prospective view 

depxctxng portions of apparatus in accordance with on! 
embodiment of the invention. 

^^^..JXL aU:in;a':~7,' 

illustrated in Pig. i. P<«tion of the apparatus 

aepicting "u^^^tr pert\o„s\/r."°""' 
35 Pios. 1 and = P^'i"""" the apparatus depicted in 



ste»s T ' " """"^ -I'Piotin, certain 

Wntl according to «.bodi»ents of the 

invention. Pig„res s and e are a diagra^S: 
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perspective view depicting portions of apparatus in 
accordance with a further embodiment of the invention. 

BE8T MODE OP CARRYING OUT IKVBMTXON 

Apparatus in accordance with one embodiment of 
5 the present invention includes three generally helical 
transmitter coils 10 disposed in a common plane. Coils 
10 are mounted in fixed position in the frame of 
reference of patient receiving bed 12. That frame of 
reference is denoted by a Cartesian coordinate system X, 

10 Y, Z as shown in Fig. 1. A patient P may be positioned 
on the patient receiving bed. The axes 11 of the coils 
are parallel to one another. A patient receiving bed 12 
extends just above the plane of coil 10. The apparatus 
further includes an object or probe 14. The probe is 

15 adapted to be inserted into a medical instrument such as 
a catheter 16 and positioned at a desired location in 
the catheter, such as at the distal tip of the catheter 
or at another location along the length of the catheter. 
Probe 14 has mounted thereon a sensor 18. Sensor 18 

20 includes three component sensors 20, 22 and 24 adapted 
to sense components of the magnetic field in mutually 
orthogonal local directions X», Y« and That is, 

component sensor 20 is sensitive to magnetic fields 
directed in direction X« , but largely insensitive to 

25 fields in directions Y« and Z*, whereas sensor 22 is 
sensitive only to fields in direction Y» and sensor 24 
is sensitive to fields in direction Z». These sensors 
are adapted to provide separate sensor signals 
representing the separate components. Sensor 18 may be 

30 a solid state sensor of the type described in the 
aforementioned International Patent Publication 
WO95/09562, the disclosure of which is hereby 
incorporated by reference herein. As fvirther described 
therein, each of the component sensors may include a 

35 generally planar magnetically sensitive film, such as a 
magneto-resistive film or a Hall effect sensing film. 
Each such film may be sensitive to fields directed 
normal to the plane of the film. Alternatively, sensor 
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18 may include an array of miniature coils, the axes- of 
the coils being oriented orthogonal to one another. 
es!ent!\, "P^-ent the preferred sensors, 

essentially any other magnetically sensitive device may 
5 be employed as, for example, magneto-optical sensors and 
rlux gate magnetometers. 

component sensors 20, 22 and 24 are connected 
through a cable 26, with separate leads for each sensor, 
to a command unit 28. command unit 28 (Fig. 3) includes 
10 an input amplification and analog to digital (..a/D") 
conversion section 30 adapted to receive the individual 
signals from component sensors 20, 22 and 24 of sensor 
18, amplify the same and convert the same to digital 
form. The amplification and a/d conversion unit 30 may 
also xnclude other conventional signal processing 
devxces such as analog or digital band pass filtering 
and noxse rejection devices and signal averagers, 
command unit 28 further includes a computation unit 32. 

Computation unit 32 mav be inni omi=>r.-t-»^ 

20 implemented as a programmed 

20 general purpose computer. As further discussed below 
the position computation unit is arranged to compute the 
dxsposxtxon Of sensor 18, and hence the disposition of 
the object or probe 14 at the catheter tip from the 
sensor signals. As used in this disclosure, the term 

25 "dxsposxtxon" of an element refers to the position of 
«xe element, the orientation of the element or both. 
Thus the computation unit is arranged to calculate the 
position of sensor 18, the orientation of the sensor, 

30 unit^'f '''k' "^r^ orientation. Command 

for ' " ^'^^'^^ -^o-) 

for provxdxng a human intelligible representation of the 

xntellxgxble representation may be provided either as 
numerxcal information presenting the position and/or 
35 orxentatxon of object 14 in the x, y, , coordinate 
system or, preferably, as a pictorial representation of 
the object and of the associated catheter superposed on 
a pictorial representation of the patient 
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Comaand unit 28 further includes a control 
unit 34. Control unit 34 is linked by output lines 36, 
38 and 40 to three separate coil drivers 42, 44 and 46. 
Each coil driver is linked to a separate one of 
5 transmitter coils 10. Each coil driver is adapted to 
send a direct current through the associated transmit 
coil 10. Each coil driver is arranged to control the 
amplitude of such current, and to turn the current on or 
off, in response to control signals received from 
10 control unit 34. The control unit is arranged to signal 
the coil drivers to provide currents to their respective 
transmit coils in alternating sec[uence, so that coil lOa 
receives current while coils lOb and 10c are inactive; 
coil 10b receives current while coils lOa and lOc are 
inactive and coil 10c receives current while coils 10a 
and 10b are inactive. The control unit receives data 
from amplification and conversion section 30 and, as 
further discussed below, actuates the coil drivers to 
vary the amplitude of the current to each coil. The 
control unit may include conventional interface devices 
such as digital to analog converters or bus interface 
unitss so that the output of the control unit is 
compatible with the control input of each coil driver. 
Also, although the control unit is illustrated 
25 separately from the other logical units of command unit 
28, it should be appreciated that the control unit may 
share physical elements of the command unit and other 
elements. For example, where the command unit 

incorporates a general purpose computer, the processor 
of the computer may serve both as an element of the 
position computation unit and as an element of the 
control unit, executing functions appropriate to the 
different units at different times. 

In a method according to one aspect of the 
invention, an catheter 16 is advanced into the body of a 
patient P. The probe 14 with sensor 18 thereon is 
disposed at the tip of the catheter. The catheter tip 
is disposed at an unknown location somewhere above the 
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^tvL"'/""' 34 actuates the coll 

drxv«e in sequence, using an initial or default value 

xo. *™plifioatxon and conversion unit 30 samples the 
s alcmal ,roB each of the co^onent sensors 20. 22 and^I 

flow tr^'i''^'" "~ '^'^^ ™^ 

flow through each coll. por e«u.ple, at a preselected 
tine after the beginning of current flL t^^ 
transmit coil loa, unit 30 taKes a sample of the s±^l 
.0 from each of component sensors 20, 22 and 24 „d 
converts the same to digital format. comrn^ JL'it 28 
then calculates a total field magnitude based upoTte^ 
lnd.v.aual Signals. The total field magnitude is: 

" |fi..»l = 

where: Bio. Is the magnitude of the magnetic 

^cS^::: " ™ " " ^- 

stren«„ " " "^^^^^^'J' ^^'''^''r relating the signal 

aurmg thettu^tiot; J^"^' ™ " 

" ^ » ^ ^""^ ''^ ^" sensitivity 

constants and signal strengths for the other sensors 22 



30 



3S 



and lo» • ^ «y=tem actuates colls lOb 

H^e aLln" °L ' "f"' ™* ^^--'th. 

t^l f magnitude of the 

total field vector prevailing at the sensor during 
actuation of coil lOb and independently computes^ 
-a^xtude Of the total field vector prelaill^ at ^e 
sensor during actuation of coll loc. 

AS Shown In Pig. 4, after detecting the 
component signals and calculating the tot^l fl^^I 
~^xtude prevailing at the sensor during actuation ^f 
each coil, the control unit determines whether or not 
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all Of the field magnitudes are within a preselected 
range of magnitudes* This preselected range is chosen 
to lie within the optimum operating range of sensor 18. 
Thus, the minimum field magnitude is selected to lie 
5 well above the noise threshold of the system and above 
the minimxim sensitivity level of the sensor, where the 
maximum field level is selected to lie well below the 
maximum limit of linearity of the sensor, and well below 
the maximum field which the sensor can tolerate without 

10 loss of accuracy. For a typical magneto-resistive 
sensor, which is most accurate and repeatable when used 
with fields of less than about 4 Gauss, the preselected 
range of field magnitudes may be from about 1.0 to about 
2.5 Gauss. With typical Hall effect sensors, which are 

15 most accurate when used with fields over about 30 Gauss, 
the preselected range will be above about 30 Gauss. If 
the three field vector magnitudes observed during 
actuation of the three units lOa, lOb and lOc all lie 
within the preselected range, than the system proceeds 

20 to calculate the position and orientation of sensor 18, 
and hence the position and orientation of probe 14 and 
the associated tip 16 of the catheter using conventional 
position finding algorithms. For example, the 

mathematical methods disclosed in United States Patent 

25 4,710,708 for finding positions using a multiple 
transmitting or receiving stations and a multi-axis 
sensor can be employed. The disclosure of said ^708 
patent is hereby incorporated by reference herein. 
Briefly, the magnitude of the fields in each of the 

30 local or sensor directions X\ Y», represented by 

each of the component sensor signals from each of the 
component sensors 20, 22 and 24 is a function of the 
overall strength of the field from the coil (also 
referred to as the magnetic dipole moment of the coil) , 

35 the distance from the particular coil to the sensor and 
the sensor rotation angles, i.e., the angles between the 
local directions X', ¥• and Z' and the coil frame of 
reference directions X, Y and Z. When three component 
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a.m.or readings are collected during actuation ot three 
separate coils .»J equated with the component strength 
jessed as a function ot location tor fields frot^; 

Iir«r t. " °^ «-= "nsor 

and the three rotation angles) . The derivation of these 
e^ations is set forth in Appendix That sy«» "f 

equations can be solved by iterative methods sul as the 
method Of Marquardt or the method of Broyden for H^l 
.0 squares solution of an overdetermined sy.t« of 
nonlinear equations. The command unit then provides 
output indicating the position and orientation of 1^e 
sensor and hence indicating the position and orientation 
Of probe 14 a^J of the distal tip of the catheter. 

" °' *^ -aanitudes is out 

Of the preselected range, the system does not calculate 
the position and orientation. Instead, control ™,it ,4 
changes the field strength of the coll or coil. 



15 



^ associated with the out-of-ran,e field strength. For 
» exa^le, where probe 14 and sensor « are r^ative^ 

" r^.^L.*"! .^-- detected 



ZZ. r ^^'"^ ^""'^ preselected 

Z7Z 4=T ' " ""^ -"-and ooil 

driver 42 to reduce the current to coU loa on the ne=ct 

c^cl- conversely, if the probe and sens^ 

sens^ d 

sensed during actuation of coil loc with the default of 
current value win be below the preseiected ran^e! 
control unit 34 therefore will instruct coil driver « 
30 to increase the current to coll loc on the „e" 
actuation cycle. The coil ^^i„. 

... , . drivers may be arranged to 

vary the overall field strength or dipole moment of each 
coil stepwise, as by varying the current stepwise. Each 
increase or decrease commanded by control unit 34 may be 
one step. Alternatively, the co«:rol unit can calculate 
an increase or decrease proportional to the degree to 
Which the total field magnitude deviates fro. aT^« 
value Within the preselected range. Thus, a relat!^" 
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large change can be made when the field magnitude is far 
outside of the range, whereas a smaller change can be 
employed when the field magnitude is close to the range 
or within the range. The correction process continues 
5 until all field strengths are within the preselected 
range, whereupon the system calculates position and 
orientation. After the system finds values of coil 
currents which produce field magnitudes within the 
range, subsequent actuation cycles use these currents. 

10 During operation, as the physician utilizes catheter 16, 
the position of the catheter tip, and hence the position 
of probe 14 and sensor 18 may change. Any such change 
may bring one or more of the field magnitudes outside of 
the preselected range, whereupon the system will 

15 readjust the currents to the coils once again. 

As the feedback control unit readjusts the 
currents to the coils, the altered values of current are 
translated into new values for field strengths from the 
individual coils which are used in the aforementioned 

20 position determining equations. In this manner, the 
system assures that the sensor is always exposed to a 
field having a magnitude within the preselected range 
wherever the sensor is placed within a sensing volxime 50 
extending over a preselected region above the plane of 

25 coils 10. The exact size of sensing volume 50 will 
depend upon the breadth of the preselected field 
magnitude range and the dynamic range of coil drivers 
42, 44 and 46, i.e., the degree to which the coil 
drivers can vary the currents. The size of sensing 

30 volume 50 within which the sensor will receive fields 
within the preselected field magnitude range from all 
coils will also depend upon the positions of the coils. 
However, for a typical system having three coils spaced 
at vertices of an equilateral triangle with sides about 

35 40 cm long, the sensing volume includes a region 
extending upwardly about 60 cm from the plane of the 
coils. At the plane of the coils, the sensing volume 
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extends about 20 cm beyond the equilateral triangle 
bounded by the colls. "langie 

m the description set forth above, only one 

5 IZZ ^""'^ *° ^"'"^^ ^-"i- an^ 

location Of only one object. However, plural objects 

and Plural sensors Bay be utilized. As described i^ t^e 
aforementioned International Publication HO9S/09562 a 
fiau:l r °\ ^ utili.^.''"^c^ 

.0 and a sensor 54 of the same type as sensor le releaUbly 
connected to the sensor body. The amplification ^d 
converson unit 3o is connected to each of th«e 
additional sensors in the same manner as to sensor 18 
The tag of each sensor may be affixed to the body of i^l 

15 patient, and an image of the n.n— » ^ 

„_,.„ . . ^ " patient may be acquired 

using any imaging method which will show the desired 
structures of the patient- s body and which wlll"^! 
Show the tags 53 as, for example, magnetic reson^e 
i«gin, („Ri,, x-ray, OAT scanning or the li,ce. Zl, 

d«. T."" facilitate imaging in thl 

desired imaging method. Where x-ray imaging is used, 
tte tag may be radioopaque; where «RI imaging is used 
the tag may Include a material which has magneti; 

« wtrr distinguishable from that of 

bodily tissue. ordinarily, the image is acquire, prior 
to use Of the catheter and probe discussed ^bove. Z 
sensors 54 typically are not present in the imagl^" 
step. After the imaging step, the sensor bodies 54 an! 
sensors of the fiducial markers are attacher tV :he 
.-dy a* the positions Of the tags. Whent^! 
catheter is employed, sensor signals from sensor 18 on 

Slim "'^'^^ concurrently with 

Signals from the sensor 52 on the fiducial marXers! The 

35 ^f"^"';"^-^'- l°-tio„ and orientation of sen^: 

filld , °* " -a^netic 

fields from coils lo in the same manner as the Lt^ 

p"rir T -i«ntation of sensor TZ 

probe 14. The acquired position and orientation of the 
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fiducial markers 52 can be used to register the 
previously acquired image of the patient with the 
position and orientation data for sensor 18 and probe 
14. For example, where the position and orientation 
5 data for the sensor 18 and probe 14 are displayed as a 
pictorial image of the catheter tip 16 on a screen, the 
position and orientation of fiducial markers, as 
determined by the magnetic field locating method can 
also be displayed. The previously acquired image data 

10 may be displayed as a pictorial representation of the 
patient's structures and that representation will also 
include a picture of the tags 53. The previously 
acc[uired images of the tags are brought into 
registration with the images of. the fiducial markers 

15 derived from magnetic location by transforming one or 
the other image until the picture of the tags in the 
previously acquired image overlies the representation of 
the fiducial marker in the magnetically acquired image. 
Such registration can be achieved by manually adjusting 

20 inputs to an image display system while visually 
matching the displayed images, as described in the 
aforementioned "562 International Publication, or else 
by automatic calculation of the transformation and 
rotation parameters necessary to bring the images into 

25 registration. Such parameters can be calculated by 
equating the positions of three or more fiducial marker 
sensors in the magnetic location frame of reference (the 
X,Y,Z coordinates of Fig. i) with the positions of the 
tags 53 in the image frame of reference as modified by 

30 a matrix of unknown rotations a further matrix of 
unknown transformations. As further described in 
Appendix B, the resulting matrix equation yields an 
overdetermined system of nonlinear equations which 
include the rotation angles and transformation 

35 distances. When the images are properly registered, the 
position of the catheter tip is displayed in true 
location and orientation relative to the patient's 
internal structures. In similar fashion, a greater 
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n^bar of fiauotal markers, or „ore than on. nadloal 
instrumant may be located. "naaical 

" "^"^ multiple sensors, such as 

sensors 18 and 54, the coils may he operated i; separa" 
cycles so as to provide each sensor with fields havin! 
rrr'V" """^-'"^ .ppropriate%or 

::rLver,rr:;d;r."dr " -^^ ~ ^ 

lo ^ ' ^ i^^noe transmitting colls 

first , '""^ "^""^ -^ing 

riiris- f-: ^ tj-^^^ 
r trzr%:iir m~ r\™- 

«ithin the preselected range a. drsc^:!:: ahov:.^"^:!^;! 
the second cycle, the signals from the sensor 54a^f ! 
first fiducial marker are acquired, the signals fro^ 
sensor 13 and the other fiducial marker sensorrs/h, 

^"virTfieT: adjusts; 

ua^iixi^uae at sensor 54a within ^-i-e 

sensor 54b are acguired whereas the other signals are 
xgnore, in the fourth cycle, the signals from srsor "c 
are acquired and others are ignor«l. The cycles Zy ll 
secruential, with the entirety of a first Cle^i^! 
performed before the entirety of . second c^cle. Z 

ZZl^' T: ^ -> anoLer. Zl 

sample, =oil driver 42 can be actuated first to provide 
a ™t appropriate to yield the proper field ^^^^ 
first s«»or 18 and then can be actuated to provide th! 

s!bT„: 54c"T """"'^^ "-^"^^ slnsors %!' 

ana 54 c in snor^Ac^^^ ^ ' 



54b and 54C in succession, followed by a simila; 
Multiple aotuation of the other coils. .ikewise , r^e" 
more sensors are employed, there may be a. many separate 
» ofihe I " "tematlvely, if 

r. t.!':"!"^^" - -ther. or 



^^4. ^ — to one another, or 

less than the number of sensors. Thus, the currLts 
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used on a given cycle may be adjusted to provide fields 
within the desired range of magnitude for one sensor* 
Signals from other sensor can be acquired during this 
cycle, although such other sensor may encounter a field 
5 outside the preselected range. In some cases, the coil 
currents are fixed for some cycles and adjusted for 
others. For example, where numerous fiducial markers 
are employed, it may be possible to achieve acceptable 
registration accuracy using fixed coil currents to 
10 provide the fields when data is acquired from the 
fiducial markers, whereas the coil currents may be 
adjusted in the manner discussed above when data is 
acquired from the sensor in the active device or 
catheter • 

15 Certain sensors tend to lose accuracy when 

exposed to magnetic fields above a predetermined 
maximum. For example, certain magneto-resistive sensors 
temporarily lose accuracy if exposed to magnetic fields 
above about 4 Gauss. Where such sensors are employed in 

20 a multi-sensor system, the system should have 
appropriate provisions to avoid exposing sensors to 
excessive fields. For example, if sensor 54a is 
disposed close to coil lOb, whereas sensor 18 is 
disposed remote from coil 10b, sensor 54a may be exposed 

25 to extremely high field levels if the system adjusts the 
current in coil 10b to the level required to generate 
the appropriate field at sensor 18. To avoid this 
condition, the data from all sensors can be acquired 
during all of the cycles, and control unit 34 can be 

30 arranged to increase coil currents in a progressive 
manner over several cycles when an increase is required. 
The data from sensors which are not used for position 
monitoring on a particular cycle can be used to inhibit 
further increases if the field at the nominally unused 

35 sensor is approaching dangerous levels. Thus, if the 
system is in the process of progressively increasing the 
coil current in coil 10b to provide an adequate field 
level at sensor 18, the system may terminate such 
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increases if field «>agnitude at second sensor 54a 

during the reading cycle associated with first sensor ^8 
reaches the maximum level allowed at i->,o ! 
If this f.nr,rt<4.- aiJ-owed at the second sensor. 

"iJ-s conditxon occurs whiio -t-H/^ -i^ , 
« ^^ui6. wnixe the field level at fir-*!*- 

, "»ors use a biasing magnetic field 

If exposed to excessive fields 

« ^ ixexas, such sensors oan k^^ 

reset and restored to accuracy after 

within the sensor. The command unit k «PP-«-iea 

trifT«o*- "«iana unit may be arranged to 

trigger the reset process far- ««« y « to 

*^ tor one sensor ^■f n- 

in the foregoing discussion, the coils ar« 
driven alternately in a ti««. ™»-.4- - , 

n^" ---- r ^^r 

Thus, each coil can be actuated with a driving current 

-nso. „ui incxuae component. vViT^t 

««ere„t carrier £.e,ue„ci.., which cnle ^ 
^ analog or aigital filtering methods. *S 
signals fro» component sensor 20 22 and / 
■0 th. oarrier freguency of coil lo. "1 k 

calculate «.e total fiela ^^it.ara\t:::L,rie T oon 
10a. That value in turn can be used to i-^i / 
control unit 34 to adjust the c^^^t to coTlo'a^ 

^> :™si" te°"" ^''"^ - .a,ustea\rsi„i:^ 

carriL-^Tre^Lrr""" ^."^ ,t 

s«,sors. Por example, each coil can ^ .^"'^ 
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multiple carrier frequencies, equal to the number of 
sensors, the carrier frequencies for each coil being 
different from one another and different from all of the 
other carrier frequencies for all of the other coils. 
The components in the signals from the first sensor 
varying at the respective first carrier frequencies can 
be detected and the strength of the coil currents at the 
respective first carrier frequencies can be adjusted 
accordingly, so that the field magnitudes at the first 
carrier frequencies prevailing at sensor 18 will be 
within the desired range. In this stage of the process, 
the components varying at the other carrier frequencies 
are ignored. The process is reversed with respect the 
fiducial marker sensor 54a. This approach can be 
extended to greater numbers of sensors. Also, the 
frequency multiplexing and time multiplexing approaches 
can be combined. Thus, each coil can be driven at only 
one carrier frequency • During a first cycle, first 
sensor 18 is actuated and the components varying at all 
of the carrier frequencies are monitored. The coil 
currents are adjusted accordingly to bring the fields at 
each carrier frequency into the desired range. During a 
second cycle, the second sensor 54a (the first fiducial 
marker sensor) is used and the coil currents are 
adjusted in the same manner, and so on with the other 
fiducial marker sensors. 

In a further variant, the coil currents, and 
hence the strengths of the fields from the individual 
transmitters, can be adjusted to bring the field 
components detected by each individual component sensor 
to within a preselected range. in such a system, each 
coil is adjusted separately with respect to each 
component sensor. In a first cycle, the current to 
transmit coil 10a is adjusted to bring the individual 
sensor signal from component sensor 20, representing the 
field component magnitude in the local X» direction, 
into a preselected range. During this cycle, the 
signals from the other component sensors 22 and 24 are 
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disregarded. In the next cycle, coil loa is readjusted 
to brxng the fieia component in the „iae pri»e direct!!^ 

5 ran" I " ~-Pon.nt sensor 22 into the optimua 
range. The sane coil loa then brings the field 
component in local direction z- into I pres.lect^a 
range, and the signal fro„ the pri„e directio^ oJ^Tr: 
section 24 IS -onitored. This sec,uence of op^rationHs 
then repeated again for earsh „t .k perations is 

.0 coils. Xhis apploach cln Tj:^:^! ttTlLaTrr 
each having plural component sensors. HetT agaiT 
systen keeps tree of the coil current »agnlt!des' u~d 
to produce the fields since bv .»,->, '"^^"des »sed 

Th. _., ^ . «>y each component sensor. 

IZZL °' ^ si^ltaneous 

Where the local direction associated with a 
particular component sensor is orthogonal, or Crly 
orthogonal, to the direction of the field produced bj a 

b^I^^tV " "'^ " ^ Lpossltie^^ 

bring the component in that local direction into the 
preselected rang, without either exceeding the curr^t 
capacity of the coil ariver or producing a totai™a 

2S Z: Z °~ ^= sensors. x" 

instance, however, at least 



30 



35 



least one of the othor- 
component sensors will receive a component hlvi^ I 
magnitude in the preselected range of magnitudej 
this variant, the signal from all of th. component 
sensors may be monitored during a cycle associated 

Lin r r 

other unused component sensor to an excessive field 
component in its sensing direction. 

m a further variant of this approach 
preselected range of magnitudes for th. f i.7d^^n^: 

a single preselected value, preferably within th. 
optimum range of accuracy of the particular ?ol;on!::: 
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sensor. The feedback control system thus operates to 
adjust the coil currents until the field component 
magnitude is at such single value. The position and 
orientation are calculated in the same way as discussed 
above- This variant has the advantage that nonlinearity 
in the component sensor response cannot affect the 
accuracy of the system. Provided that it is known that 
a particular reading from the sensor corresponds to the 
preselected value of field component magnitude, 
deviation from a linear relationship between component 
magnitude and sensor reading at other values of the 
component magnitude will not impair the accuracy of the 
system. 

Certain field sensors exhibit so-called "off- 
axis sensitivity". That is, the transfer function or 
relationship between field component magnitude along the 
sensitive axis of a particular component sensor and the 
reading from that component sensor varies when a strong 
field component orthogonal to such axis is present. 
Off-axis sensitivity can be corrected by using the 
readings from two component sensors to evaluate the 
magnitude of the field perpendicular to the sensitive 
axis of the third component sensor, and using that 
magnitude to determine a correction factor to be applied 
to the reading from the third component sensor. 

In the system described above, when the coils 
are driven with direct current, amplification and 
conversion unit 30 samples the data from each component 
sensor at a preselected delay time after the inception 
of the current. In a variant of this system, the 
command unit adjusts the delay time from inception of 
the current to sampling depending upon the current 
applied in a particular actuation. After inception of 
the current, the current gradually rises to its steady 
state value, and hence the overall field strength or 
dipole moment likewise rises gradually. The changing 
field induces eddy currents in electrically conductive 
materials which may be present in or near the sensing 
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volum.. These currents mauce additional .agnetic 
fxelds. Which can induce erroneous sensor readings ^^e 
ed^ currents decay when the rate o. change in the ii^^d 
slows as the field approaches its steady state value 

TJT ^"""^ """"^ ---- 1; 

rise to an appreciable portion of its steady state 
value, and long enough for eddy currents to dissipate 
sufficiently that they do not cause apprectahle e^ors 

'^'-^ «»xxer'"nid 

The system illustrated schematically in Figs 
5 and 6 IS Similar to that discussed above, except^"; 

^^:t" s" tf f"-^"- and sensors are ZJT. 
That IS, the probe or object 114 to be tracked is 
» equipped With a transmitters lis incorporating Lee 
miniature coils 120, .22 and X24 disposed on It^ii^ 

and coil drivers (not shown) similar to those diseased 
above with reference to Pig. 3 Th. ft^ll 
» .r»e system has three undir\cti;nars'en ^ 11"'"^:: 
and HOC mounted in fixed position relative to "e 

Cror^'"'""^ censors are mounti iH 

common plane. The common nlane i-k- 

Piane of the sensors extends 
generally vertically on one side of ^-k ^""^"^^ 

" ----- - sensing volume-Ho :Lenrou^^t 

r:::ivi'^ IT" "'^"^ " ^ --^an^ 

This system can be used in essentially the 
same manner as the system discussed above. Here aga^^ 

«<»«.1: of the transmitted field associated with elch 
coil is adjusted so that th. magnitude of I^e f Hid 

~"g a"t - --a" 

« or:::,;!"^:. " ^""'"^ ^ I-aselect«, range 

features Tsrs^dTbtrrn Te 'I^lr"" 
.aparting from the present^^ventior"": iTZl 
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example, the transmitting coils 10 of Fig. 1 can be 
disposed in a vertical plane arrangement as illustrated 
for the sensors in Fig. 5. Also, the number of 
transmitters and sensors can be varied. For example, as 
disclosed in the aforementioned International 
Publication WO94/04938, the system may include a sensor 
having three mutually orthogonal sensing axes and a set 
of three mutually orthogonal transmitting coils having a 
common center in the fixed frame of reference. 
Alternatively, the system may include a single uniaxial 
transmitting coil on the object to be located and three 
sets of sensors on the affixed frame of reference, each 
such sensor set incorporating three receiving coils or 
component sensors adapted to detect fields in a mutually 
orthogonal direction. The reciprocal arrangement — 
with three sets of mutually orthogonal transmitting 
coils and a single uniaxial sensor on the object to be 
tracked — is also usable, in general, the sensors and 
transmitters should define a plurality of transmitter- 
receiver pairs, each including one element on the object 
and one element on the fixed frame of reference. 

In a further variant, the system adjusts the 
output of the transmitters in response to the calculated 
disposition of the object being tracked, rather than 
directly in response to the component signals or total 
field magnitude signal. Thus, the system can initially 
operate with the default of current values; derive an 
initial reading of the position and orientation of the 
object and then use that initially determined position 
and orientation to calculate the desired setting for 
each coil to achieve the desired field levels at the 
sensor. Such desired setting is selected to yield the 
field within the desired magnitude range at the object, 
asstming that the object has the position and 
orientation found in the, initial reading. On the next 
cycle, the so-calculated coil currents are utilized and 
the process is repeated. In a variant of this approach, 
the system can store a lookup table listing the 
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object pos.txon .„a orientation. osing the initially 
deter„„ed position ana orientation! the s"t^ 
retrieves appropriate coil-c=urre„t ^luj trJ T 
» I-HUP table for use in the next cycle 

is as3o=i«ea"4ranrr"Te* sa""^^' 

-a With other .eaical i^stru^rn s "IT^^T 

enaoscopes ana surgical instruments i*e L.^^ 

other than neaical instruments. Por examnl= ""^^^ 
usea in tracHin, an input aevice f or a c^tlr ^ 

Of the f eaTur«"aesrib:?"h combinations 
15 ^ ^. ^^^""^^ described above can be utilized without 
15 departxng fro. the invention, the foregoing descri^i:^ 
Of the preferred embodiinents should be takL 
illustration rather th«« k ^ "^^^ °^ 

inv««.H4 ^ limitation of the 

invention as defined by the claims. 

IMDnHTpyiiT. APPI.TffiinTT -TTY 

20 

The invention can be applied in the medical 
and surgical fields as well as 4« ^-v J "^'^^^^^ 
industry. care 
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Appendix A: calculation of Position and Orientation 

5 Provided that we know the physical 

configuration of the field generator (transmitter) 
which is fixed in position during operation, the 
magnetic field detected by a sensor is a function of the 
position and orientation of the sensor. In our system 
10 field generator coils are stimulated sequentially. The 
field sensed by a sensor (3 component sensors per probe) 
can be expressed in terms of position in x, y, z and 
orientation a, y (roll, pitch and yaw respectively) 

i.e. : 

15 

B[sensor] [coil] « f [sensor] [coil] (x, y, 

Z/ OL, y) 

Where [sensor] designates a particular sensor 
and [coil] designates a particular transmitter coil. 
20 If the real field the sensor measures when 

[coil] is on is B« [sensor] [coil] , then theoretically, 

B» [sensor] [coil] = B[ sensor] [coil] 

25 

i.e. 

B' [sensor) [coil] - f [sensor] [coil] (x, y, 
2/ P, y) = 0.0 

30 

since we have 3 sensors and 3 coils, the total 
equations are 9 with 6 unlcnowns (x, y, z for probe space 
location, a, P, y for its orientation) . By applying non- 
linear least square method, we can find unique x, y, z, 
35 a, P, y for the probe. 

Above shows the general ideal of the 
algorithm. In detail: 
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Assume the orthogonal x,Y,z reference 
coordinate syste. (magnetic location Cartlsiln 
coordinate) is described by matrix Cartesian 



^112 em 
6ni em em 
.em em em. 



the probe's orthogonal system is: 



10 



Qpw Cpn epii 

Sp2i epQ e^a 
epix ep3i ep3i 



IS 



Ho ^''''^^ sensors on probe may not 

be orthogonal to each other thoi,. 

can described as: ' ^^"^ non-orthogonal axes 



^"11 e«i2 e„i3 

fi»>21 e«22 en23 
.C«3I gn32 ^1133 



20 m 1 ^ ^ transfer matrix T[i][j] which is to be used 
20 m later on calculation can be obtained from: 

Tviin~ e,ni]*efifi 

Vi.7e{U^} 

used.l« ortho_0V[i]cj] which is to be 

used also can be defined as: 



ortho OV =p 
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since we use roll (a), pitch (P), yaw(y) to 
define probe orientation, ortho_OVti] [ j ] can be also 
described by: 



10 



20 



25 



sinO) sin(y) 
sin(P) cos(Y) 



15 sinO) sin(Y) 
8in(P) cos(Y) 



ortho_OV[l] [1] = cos (a) cos(y) - sin (a) 

ortho_OV[l) [2] = cos(a) sin(Y) - sin(a) 

ortho_OV[l] [3] = -sin(a) cos(p) 

ortho_ov[2] [1] = -cos(p) sin(Y) 

ortho_OV[2] [2] « cos(P) cos(Y) 

ortho_OV[2] [3] = sin(3) 

ortho_0V[3] [1] = sin (a) cos(y) + cos (a) 

ortho_ov[3] [2] = sin (a) sin(Y) - cos (a) 

ortho_0V[3] [3] = cos (a) cos(P) 



orthogonal vector matrix can therefore be 
calculated by matrix multiplication of previous defined 
matrix T and orhto__OV: 

. .. ov-T*artho_OV 

The theoretical magnetic field for an 
orthogonal system generated at sensor position pointing 

to ^I'l direction can be expressed as: 

30 /lcoir\[i'\(x,y,z,a,fi,r) 

(detail abbreviated) , where jf is a known function, and 
includes a dipole moment term having magnitude 
proportional to the current flow in the particular coil. 
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. non_orthogonality correction (sensors 

not be perpendicular to each other,, the magnetic 

field sensor measures should be- y 

5 



Assume the real field sensor detected when 
10 coil is on is B'isensor][coil\ 

I unen : 

a'tKBSorHoM/] - S[«mo7-]ti»(/] » 0.0 

therefore, the 9 «,uatlo„s to be solv«J for , 
« y, z, a, p, 1, arei 



25 



S'[«««»-2]I<^/>J-|/t<:<,,n][,Kr.j-.z.a.Ar)xoV-'[0[»™<-2) = 0.0 
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3 

=0,0 



5'[5e«sor2][co//3]- ]^/[co//3][i](x,>^,r,a,Ar) x ov*'P][^^'wor2] = 0.0 



1=1 

5 



B'[sensor3][coin] - ^f[com][iXx,y,z,a,M x ov'^PlE^^w^^^'-S] = 0.0 



5'[^ensor3][co//2]^ ]^/[co//2][/](x,>;,z,a,^,r) x ov"'W[^^w^or3] = 0.0 



Is] 



10 

3 

5»[5en^or3][co/73]-- ]^/[coi73][jl(x,j/,r,a, Ay) x ov''P]['y^'»-y^''3] = 0.0 



15 



A well-known non-linear least squares equation 
solver can be applied to solve the above equations and 
find probe position y, z and orientation fi, y. 
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Appenaix Bt Registration of Piduoial Markers 

Asstunption: 

. . coordinates of three fiducial 



10 My3, Mz3)} 



{(Mxl, Myl, M21), My2, M22) , (Mx3, 



15 



^'^'^ coordinates of the same three 

fiducal points in tho v v .7 ^. ^^ee 

^ * X,Y,2, coordinate system of the 

fixed frame of reference r^f 

apparatus: -^-^^nc^ of the magnetic locating 



Py3, Pz3)> 



{(PXI, Pyl, Pel), (p^2^ p^^j^ ^p^^^ 



20 



25 



(X3, y3, 23)) can be formed by { (ai, bl, ol) (a2 b-> 

s:'r:::;i"' -Vn-sxrion": 

Tne rotation matrix R is: 



iixr i?)af /?zr 
^ Rxy Rzy 
Hxz Ryz Rzz 



30 



where 
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Bxx = cos(a) X cos(y ) - sin(a) x anOff) x ^n(y) 
Rxy = cos(a) x anO^) + an(a) x sin(^ x cos(y ) 
Rxz = -an(a) x cosO?) 
Ryx = -cosO?) X sin(y) 
Ryy = cos09) x cos(y) 
Ryz = an(^ 

iizx = an(a) x cosO^) + cos(a) x anQS) x sinO^) 
= sin(a) x anO^) - cos(a) x sinO?) x cosO^) 
Rzz = cos(a) X cos(fi) 

The -translation matrix T is: 



T = 



X y 2 
X y 2 
X y 2 



The relationship between image position and 
magnetic location apparatus frame of reference position 
is: 



10 



Mxl 


Myl 


AdzV 




'Pxl 


Pyl 


Pzi 


Mx2 




Mz2 


R + T = 


Pxl 


Py2 


Pz2 


Mx3 


My3 


Mz3_ 




Px3 


Py3 


Pz3 



or 



Mxl Myl Mz\ 
Mxl My2 Mz2 
Mx3 My3 Mz3 



R + T- 



Pxl 


Pyl 


Pzl 




"o 


0 


0' 


Pxl 


Py2 


Pz2 




0 


0 


0 


Px3 


Py3 


Pz3 




0 


0 


0 
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This forms 9 single equations: 
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A/ocl X 72xx + iWyl X i?xv + A/zl X /2x2 + X - Prl = 0 
A&l X + A^l X i?^ + X i?^ +^ _ = 0 
A&l X iJzx + AO/l X ijg, + A&l X iJzz + z - = 0 
A&2 X /?xr + Afy2 X /tcv + X + X - i>x2 = 0 
Afr2 X + A^2 X Ryy + Aiz2 xRyz+y-Py2 = 0 

Mx3 X ibof + A4y3 X /fccy + Afe3 X iter + X- i>r3 = 0 
Mx3xR2x + My3xRzy + Mz3xR2z + z-Pz2 = 0 



with 6 unknowns: a. B v ^n-r- -^r.*-^*.* 
5 «^'P/Y for rotation angles, 

5 X, y, 2 for space translation. 

The registration procedure is to find rotation 
angles a, y and space translation x, . ^y solving 
the above 9 equations. ^ 

10 

in pseudo-code, the registration procedure is: 

registration () 
{ 

Initialize data buffer; 

Assign fiducial reference points to image 
data (usually 3 points); ^ 

Initialize magnetic locating system; 
Measure lab positions ^-v^ 
M) fiducial points; sxtxons of the chosen 

Apply equation solving method (non-linear 
least square) on equations (*) , 

=pace ^ansxation x,"" «' ^ - 

- Fo™ rotation „atri. « ,.0. a. p, , 

translation matrix T from x, y, z; ' P. Y and 
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The rotation matrix R and translation matrix T 
can be applied to subsec[uent positions obtained, 
including the positions of the of the sensor associated 
with the catheter. 
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CLAXliS: 

1 . Apparatus for determining position 
comprising: 

(a) field generating means (lo) for 
producing a plurality of magnetic fields, each said 
field having parameters varying with location within a 
sensing volume (50, according to a known pattern of 
variation, said patterns of variation being different 
for different ones of said fields; 

(b) at least one sensor (18) adapted to 
detect one or more of said parameters prevailing at such 
sensor when said sensor is at an unknown location within 
said sensing volume and to provide one or more sensor 
signals representing said one or more parameters; 

(c) calculation means (32) for 
determining a calculated location of said sensor based 
upon said sensor signals and said known patterns of 
variation; ana 

«3) feedback control means (34) for 
adjusting said «eld generating means to alter said 
known pattern of variation of at least one of said 
fields responsive to said sensor signals, said 
calculated location or both so as to maintain said 
parameters of each such altered field prevailing at said 
sensor within a preselected range. 

2. Apparatus for determining position 

comprising: 
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(a) a plurality of sensors (110) adapted 
to detect one or more parameters of a magnetic field 
prevailing at each such sensor and to provide one or 
more sensor signals representing said one or more 

5 parameters, said sensors being disposed adjacent a 
sensing volume (150) at different positions, different 
orientations or both; 

(b) field generating means including at 
least one magnetic field transmitter (118) movable to an 

10 unknown location within said sensing volume and adapted 
to provide at least one magnetic field in said sensing 
volume, so that the magnetic field from each transmitter 
varies in according to a known pattern of variation with 
location relative to such transmitter whereby said 

15 sensor signals will represent one or more parameters of 
the field or fields from said at least one transmitter; 

(c) calculation means (32) for 
determining a calculated location of said at least one 
transmitter based upon said sensor signals and said 

20 known patterns of variation; and 

(d) feedback control means (34) for 
adjusting said field generating means to alter said 
known pattern of variation of at least one of said 
fields responsive to said sensor signals, said 

25 calculated location or both so as to maintain said 
parameters of each such altered field prevailing at each 
said sensor within a preselected range. 
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3. 



Apparatus as claimed in claim i or claim 
2 Wherein said field generating means is operative to 
provide DC magnetic fields within said sensing volume. 

4. Apparatus as claimed in claim 3 wherein 
3 said field generating means includes a plurality of 

coils (10, 120, 122, 124) and means for intermittently 
actuating each said coil (34, 42, 44, 46, with a direct 
current, said feedback control means includes means for 
altering the magnitude of each said direct current, and 

10 Wherein said calculation means includes means for 
sampling the signals fro. each sensor when a delay time 
has elapsed after commencement of each said intermittent 
actuation, said feedback control means further including 
»eans for varying the delay time for each actuation in 

15 accordance with the magnitude of the current applied in 
such actuation so as to reduce the delay time when 

lesser cxirrents are emolov^iri 

aire employed and increase the delay time 

when greater currents are employed. 

5. Apparatus as claimed in claim 1 or claim 
20 2 Wherein said field generating means is operative to 

provide AC magnetic fields within said sensing volume. 

6. Apparatus as claimed in claim i or claim 
2 Wherein each said sensor has maximum accuracy within a 
preselected optimum operational range for each said 
parameter and said preselected range for each said 
parameter substantially corresponds to said optimum 
operational range. 

7. Apparatus as claimed in claim i or claim 

2 wherein said preselected range for each 

lor each said parameter 
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consists Of a single value, and wherein said feedback 
control means is operative to maintain each said 
parameter of the field at each sensor substantially at 
said single value. 
5 8. Apparatus as claimed in claim 1 wherein 

said field generating means includes a plurality of 
transmitters (10) disposed adjacent said sensing volume. 

9. Apparatus as claimed in claim 2 or claim 
8 wherein each said transmitter is operative to radiate 

10 a field having magnitude decreasing with r" where r is 
distance from the transmitter and n is a number greater 
than 1. 

10. Apparatus as claimed in claim 9 wherein 
each said transmitter is operative to radiate a field 

15 having magnitude substantially equal to K/r^ where K is 
a real number, and wherein said feedback control means 
is operative to adjust K. 

11. Apparatus as claimed in claim 1 further 
comprising a sensor body (18), said at least one sensor 

20 including a plurality of component sensors (20, 22, 24) 
disposed oh said sensor body, each said component sensor 
being operative to measure the magnitude of a magnetic 
field component in a preselected local direction 
relative to the sensor body. 

12. Apparatus as claimed in claim 11 wherein 
said field generating means includes a plurality of 
transmitters (10) disposed at spaced-apart locations 
adjacent said sensing volume, said feedback control 
means being operative to adjust each said transmitter to 
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transmit a plurality of fields so that for each said 
field, the magnitude of a component of such field in 
the local direction of one said component sensor will be 
within a preselected range. 

13. Apparatus as claimed in claim ii wherein 
each said component sensor is a Hall effect sensor, a 
magnetoresistive sensor, a magnetooptical sensor or a 
flux gate laagnetometer • 

14. Apparatus as claimed in claim i or claim 
10 2 Wherein each said sensor includes a coil and wherein 

said field generating means is arranged to provide AC 
fields in said sensing volume. 

15. Apparatus as claimed in claim i wherein 
said field generating means includes a plurality of 

15 transmitters disposed at spaced-apart locations on one 
side Of a reference plane and wherein said sensing 
volume extends to the other side of said reference 
plane. 

16. Apparatus as claimed in claim 15 wherein 
20 said transmitters are disposed substantially in a common 

plane. 

17. Apparatus as claimed in claim 2 wherein 
said Plurality of sensors are disposed at spaced-apart 
locations on one side of a reference plane and wherein 
said sensing volume extends to the other side of said 
reference plane. 

18. Apparatus as claimed in claim 17 wherein 
said sensors are disposed substantially in a common 
plane. 



25 
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19. Apparatus as claimed in claim 1 or claim 
2 wherein said at least one sensor includes a plurality 
of sensors, and wherein said sensors can be impaired by 
exposure to fields having components above a 

5 predetermined threshold magnitude, and wherein said 
feedback control means includes means for controlling 
said field generating means to assure that no sensor is 
exposed to a field component having magnitude above its 
threshold value. 

20. Apparatus as claimed in claim l or claim 
2 wherein said at least one sensor includes a plurality 
of sensors and wherein said sensors can be impaired by 
exposure to fields having components above a 
predetermined threshold magnitude but can be restored 

15 after such exposure upon application of a reset routine, 
said feedback control means including means for 
initiating such reset routine in the event of such 
exposure . 

21- Apparatus for sensing the disposition of 
20 an object in a frame of reference comprising: 

(a) transmitter means (10) including at 
least one transmitter for providing a field; 

(b) sensor means (18) including at least 
one sensor for sensing one or more parameters of a field 

25 prevailing at each such sensor and providing one or more 
sensor signals indicative of such parameters, said 
transmitter means and said sensor means cooperatively 
defining a plurality of transmitter-sensor pairs, each 
including one transmitter and one sensor as elements of 



wo 96/41119 

PCT/US96/08411 

-42- 

the pair, one element of each such pair being disposed 
on the Object and the other element of each such pair 
being disposed at a known disposition in said frame of 
reference; 

(c) calculation means (32) for 
determining the disposition of said object in said frame 
Of reference based upon said sensor signals; and 

(d) feedback control means (34) for 
adjusting said transmitting means responsive to said 

10 sensor signals to thereby maintain at least one of said 
sensor signals within a preselected range. 

22. A method of determining the disposition 
Of an Object in a frame of reference comprising the 
Steeps of: 



(a) 



providing one or more magnetic 

fields each having a known rtati-or-r. 

jixiown pattern of variation with 

respect to said frame of reference; 

(b) sensing one or more parameters of 
said one or more magnetic fields prevailing at said 
^0 Object to provide one or more sensor signals indicative 
of «aia one or »or« parameters; 

(e) calculating the aisposltion of saia 
Object in saia frame of reference basea upon saia sensor 
signals ana saia known pattern or patterns of variation, 



29 ana 



(d) altering said known pattern of 
variation of at least one of saia one or more magnetic 
fields responsive to said sensor signals, said 
calculated disposition or both so as to maintain saia 
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parameters of each such altered field prevailing at said 
object within a preselected range. 

23. A method of detenaining the disposition 
of an object in a frame of reference comprising the 
s steps of: 

(a) providing one or more magnetic 
fields each having a known pattern of variation in space 
with respect to object; 

(b) sensing one or more parameters of 
10 said one or more magnetic fields prevailing at one or 

more sensors having known disposition in said frame of 
reference to provide one or more sensor signals 
indicative of said one or more parameters; 

(c) calculating the disposition of said 
15 object in said frame of reference based upon said sensor 

signals and said known pattern or patterns of variation 
of said one or more magnetic fields; and 

(d) altering said known pattern of 
variation of at least one of said one or more magnetic 

10 fields responsive to said sensor signals, said 
calculated disposition or both so as to maintain said 
parameters of each such altered field prevailing at at 
least one of said sensors within a preselected range. 
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